Numerical implementation of seismic wavefield operators via path integrals
نویسندگان
چکیده
The seismic imaging problem centers around mathematical and numerical techniques to create an accurate image of the earth’s subsurface, using recorded data from geophones that capture reflected seismic waves. Using a path integral approach, a wavefield extrapolater can be expressed as a limit of depth-sliced path steps through a variable velocity medium. An image is created from the correlation between upward and downward going waves. We report on the mathematical issues that arise in implementing numerical algorithms based on the path integral approach, in particular convergence, stability, and accuracy.
منابع مشابه
Gabor Multipliers for Seismic Modeling and Wavefield Propagators
Nonstationary processes in seismic modeling may be represented numerically by Gabor multipliers and generalized frame operators. These numerical algorithms are based on a localized version of the Fourier transform and share many of the speed and accuracy benefits of the FFT, while modeling nonstationarity. We demonstrate applications to seismic deconvolution and numerical wavefield propagation,...
متن کاملConstruction of finite-frequency kernels using adjoint methods
This chapter contains excerpts from “Seismic tomography, adjoint methods, time reversal, and banana-doughnut kernels,” by Jeroen Tromp, Carl Tape, and Qinya Liu. My primary contribution to this study was to adapt a 2D SEM wave propagation code to construct finitefrequency kernels. In a series of numerical experiements, I illustrated the formation of finitefrequency sensitivity kernels via the i...
متن کاملSeismic scalar wave equation with variable coefficients modeling by a new convolutional differentiator
a r t i c l e i n f o a b s t r a c t Keywords: Wave equation with variable coefficients Heterogeneous media Accurate modeling of seismic wave Localized operator Suppression of numerical dispersion Studying seismic wavefields in the Earth's interior requires an accurate calculation of wave propagation using accurate and efficient numerical techniques. In this paper, we present an alternative me...
متن کاملA Generalized Deconvolution Approach for Local Radon Transforms
A chief problem in seismic data processing is the filtering of unwanted events like ground roll and multiples. Methods to deal with this problem often exploit moveout or curvature differences between offending events and the events one would like to preserve (primaries). In particular, removal of multiples based on moveout discrimination can be attained via parabolic and hyperbolic Radon transf...
متن کاملAdaptive finite difference for seismic wavefield modelling in acoustic media
Efficient numerical seismic wavefield modelling is a key component of modern seismic imaging techniques, such as reverse-time migration and full-waveform inversion. Finite difference methods are perhaps the most widely used numerical approach for forward modelling, and here we introduce a novel scheme for implementing finite difference by introducing a time-to-space wavelet mapping. Finite diff...
متن کامل